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ABSTRACT 

We investigate the response of the star formation efficiency (SFE) to the 
main parameters of simulations of molecular cloud formation by the colli- 
sion of warm diffuse medium (WNM) cylindrical streams, neglecting stellar 
feedback and magnetic fields. The parameters we vary are the Mach number 
of the inflow velocity of the streams, M s j n f , the rms Mach number, M Sj b g d 
of the initial background turbulence in the WNM, and the total mass con- 
tained in the colliding gas streams, Mini- Because the SFE is a function 
of time, we define two estimators for it, the "absolute" SFE, measured at 
t = 25 Myr into the simulation's evolution (SFE a b s .25), and the "relative" 
SFE, measured 5 Myr after the onset of star formation in each simulation 
(SFErd^). The latter is close to the "star formation rate per free-fall time" 
for gas at n — 100 cm~ 3 . We find that both estimators decrease with increas- 
ing A4i n f, although by no more than a factor of 2 as M s ; n f increases from 1.25 
to 3.5. Increasing levels of background turbulence (injected at scales compa- 
rable to the streams' transverse radius) similarly reduce the SFE, because 
the turbulence disrupts the coherence of the colliding streams, fragmenting 
the cloud, and producing small-scale clumps scattered through the numerical 
box, which have low SFEs. Finally, the SFE is very sensitive to the mass of the 
inflows (at roughly constant density and temperature), with SFErei^ decreas- 
ing from ~ 0.4 to ~ 0.04 as the mass in the colliding streams decreases from 
~ 2.3 x 1O 4 M to ~ 6OOM or, equivalently, the virial parameter a increases 
from ~ 0.15 to ~ 1.5. This trend is in partial agreement with the prediction 
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by iKrumholz fc McKed (|2005l ). since the latter lies within the same range as 
the observed efficiencies, but with a significantly shallower slope. We conclude 
that the observed variability of the SFE is a highly sensitive function of the 
parameters of the cloud formation process, and may be the cause of significant 
scatter in observational determinations. 
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1 INTRODUCTION 



The control of the star formation efficiency (SFE) by turbulence is a central issue in our 
present unders tanding of star formation, and currently a topic of intense study (see, e.g., 

20071 ). In recent years, several 



the reviews by 



Mac Low fc Klessen 



2004 



McKee fc Ostriker 



groups have studied the SFE of molecular clouds (MCs) using numerical simulations of 
isothermal turbulence, in which the entire numerical box represents the interior of a molecu 
lar cloud (see, e.g., the re views by 



Vazquez-Semadeni 



Mac Low fc Klessen 



2004 



Ballesteros-Paredes et al. 



20071 : 



20071 ). One confusing issue is that simulations of driven turbulence seem 



to indicate that the SFE decreases as the turbulent rms Mach number M K increases (e.g. , 



Klessen. Heitsch & Mac Low 


2000; 


Vazauez-Semadeni. Ballesteros-Paredes & Klessen 


2003: 


Vazquez-Semadeni, Kim & Ballesteros-Paredes 


2005 


), while simulations of decaying turbu- 



lence suggest that the SFE increases with increasing M s ( Nakamura fc Lil 12005 ). This has 
prompted the question of how does turbulence actually originate and behave in real MCs. To 
answer this question, it has become necessary to investigate the entire evolutionary process 
of MCs. 

The formation of MCs by collisions of warm neutral medium (WNM) streams has been 
intensely studied in recent years. A vast body of numerical simulations has shown that mod- 
erate, transonic compressions in the WNM can nonlinearly trigger a phase transition to the 



cold neutral medium 



2002 



(CNM) (e.g., 



Walder fc Folini 



Hennebelle fc Perault 



19991 : 



Koyama fc Inutsuka 



2000 . 



20001 ) . and that the dense gas produ ced by this mechanism is over 



pressured with respect to the mean WNM thermal pressure (j Vazquez-Semadeni et al. 



and turbulent, due to the combined action of Kelvin-Helmholz, thermal (iFieldl Il965l ) and 



2006|) 
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nonlinear thin-shell (IVishniac 



produced by this mechanism is continually driven for as 
The physical scenario of MC evolution was outlined by 



1994 ) instabilities (IHeitsch et alJl2005L 120061 ). The turbulence 

ong as the compression lasts. 



Hartmann. Ballesteros-Paredes fc Bergin 



(120011 ). who estimated the column densities necessary for the cloud to become self-gravitating, 



molecular, and magnetically supercritical, finding them to be compara ble. The one-dimension al 



physic al co nditions in the dense gas were c alculated analytically by 



Hennebelle fc Perault 



(119991 ) and 



Vazquez-Semadeni et al 



(200 6j). A recent review of the subject has been pre- 



sented by iHennebelle. Mac Low, fc Vazquez-Semadenil (120081 ) . 

More recently, simulations including self-gravity and "sink particles", which represent 
gravitationally collapsed objects (stars or stellar clusters), and using finite-duration com- 
pressions in the WNM, although lacking stellar feedback and magnetic fields, have been 
used to study the evolution of the tur bulent motions and of the SFE in a self-consistent 
manner (j Vazquez-Semadeni et al.ll2007l . hereafter Paper I). Concerning the velocity disper- 
sion in the clouds, these authors found that the turbulence is intermediate between driven 
and decaying, since what decays is the driving rate of the turbulence as the inflows weaken 
with time. However, they also found that the random motions are gradually replaced by 
global infall motions, as the cloud begins to contract gravitationally. Concerning the SFE, 
Paper I measured the masses of dense gas A^dense and of the collapsed stellar objects -M sta rs 
in the simulations, allowing a measurement of the SFE, defined as 

•M. stars 



SFE 



(1) 



M 

dense + M stars 

The resulting SFE was however too high, reaching ~ 50% roughly 6 Myr after the time 
at which star formation (SF) had begun (denoted tsp), although this excessive SFE can 
possibly be attributed to the neglec t of stellar feedback in that simulatio n. Indeed, Paper 



I esti mated, using a prescription by iFranco. Shore &: Tenorio-Tagld (119941 ) and a standard 
IMF ( )Kroupall200ll ). that by 3 Myr after t$p, enough massive stars would have formed as to 
be able to destroy the cloud by ionization. At that point, the SFE w as ~ 15%, closer to the 
typical values < 5% reported observationally for full MC complexes (IMyers et al.lll986l ). 

One obvious possible reason for the relatively high values of the SFE obtained in this 
type of simulations is the neglect of stellar feedback and magnetic fields. However, it is 
also of interest to perform a more systematic investigation of the degree of variability that 
the SFE could exhibit within the scenario of Paper I, simply by varying the parameters of 
the WNM compressions triggering the formation of the cloud complex. In this paper we 



4 Rosas-Guevara et al. 

undertake a first approach to such task, by varying three parameters of the WNM stream 
collisions modeled in the simulations. First, we consider the inflow speed v ; n f and the velocity 
dispersion of the background turbulence initially present in the medium, both measured by 
their respective Mach numbers, M s inf and M Sibgd , with respect to the unperturbed WNM. 
Subsequently, we consider the mass in the colliding streams -M m f, as determined by their 
radius i? in f and length / inf . Since the parameter space covered by these three parameters is 
already quite large, in this work we do not consider variations in the collision angle of the 
streams, instead having them collide head-on in all cases. 

The paper is organized as follows. In $2]we describe the numerical model and experiments. 
In §|3] we present our results, and in §H]we present a summary and our conclusions. 



2 NUMERICAL MODEL AND EXPERIMENTS 



For the numerical simulations, we use the same numerical setup as that used in Paper I, 
ex cept that we no w use the smoothed particle hydrodynamics (SPH) + iV-body code Gadget- 



2 (jSpringell 120051 ) (in Paper I we used the previous version of the code, Gadget), modified 



to include random t urbulence driving and sink particles according to the prescription of 



Jappsen et al.l (l2005f) . and including paramet erized heating and cooling, as applied in Paper 
I, using the fit of iKoyama fc Inutsukal (120021 ) to a variety of atomic and molecular cooling 
processes. This cooling function causes the gas to be thermally unstable, under the isobaric 



mode, in the density range 1 < n < 10 cm 3 . We assume that the gas is all atomic^], with a 
mean atomic weight of 1.27. The numerical box is periodic, with size L D ox- In all cases, we 



use 118 3 = 1.64 x 10 6 particles, and set the mean number of 



volume to 40. According to the criterion of 



Bate fc Burkert 



p artic 



es within a smoothing 



(119971 ) . the mass resolution 



is twice the number of particles within a smoothing volume, or ~ 80 times the mass per 
particle. The critical density for sink formation is set at 3.2 x 10 7 cm -3 , and the outer sink 
accretion radius is set at 0.04 pc. 

An initial turbulent velocity field of one-dimensional velocity dispersion, characterized 
by its rms Mach number M s t, gc j and applied at scales between 1/4 and 1/8 of the box size, is 
added to the inflow velocity field, in order to trigger the instabilities that render the cloud 
turbulent. Note that this added turbulent velocity field is applied by turning on the random 
driver for the first few timesteps of the simulation's evolution, and that what we actually 



1 Our "molecular clouds" are thus only so in the sense of density and temperature, but not of chemical composition. 
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control is the energy injection rate parameter. Thus, simulations intended to have the same 
turbulence strength do so only approximately, as the flow's response is slightly different in 
every realization. 

The initial conditions consist of a uniform medium at n — 1 cm -3 and T = 5000 K, 
in which two cylinders of length Z in f along the x direction and radius _R; n f are set to collide 
head-on at the x = Lb ox /2 plane of the simulation (refer to Fig. 1 of Paper I). Note that the 
cylindrical inflows are entirely contained within the numerical box, since the boundaries are 
periodic. The length /; n f is measured from the central collision plane, and is always shorter 
than the half-length of the box, implying that a small region between the edge of the inflows 
and the box boundaries is not given any velocity. This region is partially evacuated during 
the subsequent evolution of the simulations, as the gas within it tends to fill the void left by 
the inflows. 

Tabled] shows the various runs we performed for the present study, indicating the relevant 
parameters for each one. In addition to the inflow length, radius, and velocity, defined above, 
Table Q] gives the Mach number M Sj i n f corresponding to v ; n f at the initial temperature of the 
gas, for which the adiabatic sound speed is 7.54 km s _1 , the one- dimensional rms Mach 
number of the initial turbulent motions, M S) b g d, the total mass contained in the simulation 
box, A^box, the mass contained in the two inflows, M.- m ii arid the mass per SPH particle, 
-Mpart- The runs are labeled mnemonically, with their names giving, in that orded, the inflow 
Mach number M Sj i n f, the background Mach number M Sj b g d, and the inflow mass, .Mmf- 

Note that run Mil-Mb.02-Ma2e4 has the same parameters as run L256Aw0.17 from 
Paper I, and we take these as the "fiducial" set of parameters. Note, however, that these two 
runs are not identical because they were performed with different codes. Gadget-2 differs 
in many ways from Gadget, and in particular it takes longer timesteps, implying that the 
initial forcing (used to trigger the instabilities in the dense layer) is applied at different 
time intervals, and with different random seeds. Thus, the two runs are similar only in a 
statistical sense, but are not identical. For reference, in Fig. [1] we show a face-on view of 
Mil-Mb.02-Ma2e4 at the time when it is beginning to form stars, showing that the general 
morphology it develops is similar to that of run L256At>0.17 from Paper I (compare to Fig. 
4 of that paper, noting that the linear scales shown are different in each figure). 

We report the SFE as defined by eq. (CQ), with the "dense" gas defined as that with 
number density n ^ 100 cm~ 3 . However, because the SFE is actually a function of time, and 
sink formation begins at different times in different runs, in order to report a number for the 
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Figure 1. Face-on image in projection of run Mil-Mb.02-Ma2e4 at the time it is beginning to form stars. Note the central ring 
and the radial filaments, charactristic of this type of runs (Paper I). 



SFE, we estimate it in two different ways. One is to measure the "absolute" SFE 25 Myr 
after the start of the simulation, which we denote as SFE a b s ,25- The other is to measure the 
"relative" SFE, which we define as the SFE 5 Myr after the onset of star formation in the 
simulation, and which we denote as SFE re i 5 . This is actually clo se to the "star format i on rat e 
per free-fall time", SFRfl- (after SF has begun), as defined by iKrumholz fc McKed (120051 ). 
since the free-fall time for gas at n ~ 100 cm -3 is ~ 4.6 Myr. 



3 RESULTS 

3.1 SFE vs. inflow velocity 

In this section, we consider the dependence of the SFE on the inflow velocity of the colliding 
streams, t> m f Figure [2] shows the evolution of the dense gas mass (top panel) and the sink 
mass (bottom panel) for runs Mil-Mb.02-Ma2e4, Mi2-Mb.02-Ma2e4, and Mi3-Mb.02-Ma2e4. 
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t[U;.t] 




Figure 2. Evolution of the dense gas mass (left panel) and sink mass (right panel) for runs Mil-Mb.02-Ma2e4, Mi2-Mb.02- 
Ma2e4, and Mi3-Mb.02-Ma2e4, which differ only by the Mach number of the inflows (indicated by the "Mi#" entry in the 
run's name), having M s ; n f = 1.25, 2.5 and 3.5 respectively. 

These runs have all parameters equal, except for the speed of the inflows (cf. Tabled]), which 
are varied from M sM = 1.25 in Mil-Mb.02-Ma2e4 to M sM = 3.5 in Mi3-Mb.02-Ma2e4. Note 
that the time at which the runs begin to form sinks, tsF, is different in each case. The top 
panels of Fig. [3] then show the SFE, defined as in eq. (pQ). The top left panel shows SFE aDS) 25, 
i.e., starting from the beginning of the simulation, and up to a total time of 25 Myr. The 
top right panel shows the SFE starting from the time at which sink formation begins in each 
run, allowing one to read off SFE re i i5 . 

From the top right panel of Fig. [3] we see that the mean slope of the curve SFE(t), denoted 
(SFE), decreases, although moderately, with increasing M Sjm f. It is interesting, however, that 
all three runs form stars at roughly the same rate (i.e., the slopes of the curves in the right 
panel of Fig. [21 which give the mass accretion rate onto the sinks, and which we identify 
with the star formation rate, SFR, are all similar)!^ The decrease in (SFE) is thus due to 
the larger mass growth rate of the cloud induced by the larger inflow velocities (Fig. [21 top), 
not to a smaller SFR. 

The larger-mass clouds appear to be incapable of forming stars at a proportionally larger 



the c 
(e.g. 


ouds 1 


5 which have been shown to reduce the SFR in simulations of driven turbulence 


Klessen. Heitsch & Mac Low 


2000: 


Vazauez-Semadeni. Ballesteros-Paredes & Klessen 


2003; 


Mac Low & Klessen 


2004; 


Vazquez-Semadeni. Kim & Ballesteros-Paredes 


2005 


). This 



is illustrated in the bottom panel of Fig. [3] which shows the evolution of the one-dimensional 



2 Note that, because the cloud's mass is in general not constant in time, eq. JTJ implies that A^donse x (SFE)t^ SFR. 

3 Note that this velocity dispersion, which is the response of the flow in its dense regions to the collision of the inflows, is 

different from the initial background turbulent velocity dispersion applied to the runs, measured by the parameter M s ,bgd- 
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Figure 3. Top panels: Evolution of the "absolute" SFE (SFE a t, s .25, left), shown out to 25 Myr after the start of the runs, and 
the "relative" SFE (right), shown from the onset of sink formation, for runs Mil-Mb.02-Ma2e4, Mi2-Mb.02-Ma2e4, and Mi3- 
Mb.02-Ma2e4. SFE re i g is the value of this curve at a relative time of 5 Myr. Bottom panel: Evolution of the one-dimensional 
velocity dispersion perpendicular to the direction of the inflows for the same three runs. 

w rms perpendicular to the direction of the inflows for the same three runsfl The trends of 
SFE a bs,25, SFE rc i )5 , and v ITDS with the inflow speed are summarized in Fig. HI In this figure, 
f rms is measured for Run Mil-Mb.02-Ma2e4 at the time when it exhibits a minimum in the 
bottom panel of Fig. [3], i.e., t — 8 Myr. This is equivalent to ~ 2/3 t cros , where t cr0 s = Anf/'Wmf 
is the inflow crossing time, the time taken by the tail of the inflowing cylinders to reach 
the collision plane. We do this in an attempt to capture the true velocity dispersion in the 
cloud after the initial transients have ended, but before global collapse sets in. The latter is 
indicated by the smooth rise in t> rms for this run during the time interval 9 < t < 23 Myr. 
Then, for consistency, we measure v rms for the other two runs also at 2/3 of their own i cros . 



4 Note that the velocity dispersions in this figure arc significantly smaller than those reported in Figs. 5 (bottom panel) and 9 
of Paper I. In that paper, those figures suffered from a typographical error, having the velocities erroneously multiplied by one 
too many factors of the velocity unit, vo = 7.362 km s _1 . Thus, the correct velocities in those figures are obtained by dividing 
by this factor. 



Star Formation Efficiency in Simulations of Molecular Cloud Formation 9 



0.1 



0.6 



^ 0.4 
in 



0.2 



0.0 





- 






V -O 
rms . .- - 


- 




SFE abs g 5 






^ "A 


0'' 






- Mil-Mb.02-Ma2e4 


Mi2-Mb.02-Ma2e^ 








MI3-Mb.02-Ma2e4 



3 



> 



M 



s.inf 



Figure 4. Dependence of SFE abSj2 5 and SFE roli5 on the inflow speed for runs Mil-Mb.02-Ma2e4, Mi2-Mb.02-Ma2e4, and 
Mi3-Mb.02-Ma2e4. Both SFEs are seen to decrease with increasing i>i n f, apparently due to an increase of the turbulent velocity 
dispersion « rm s as the inflow speed increases, which in turn reduces the SFR. 



3.2 SFE vs. background turbulence strength 

We now consider the response of the SFE to the amplitude of the initial turbulent velocity 
field. Note that this field was not originally intended to produce density condensations on its 
own, but just to sufficiently disorganize the inflow velocity field as to trigger the instabilities 
that render the cloud turbulent. However, in the cases of stronger turbulence, we do observe 
clump formation everywhere in the box as a result of the initial background turbulence, and 
not just at the collision site of the inflows. 

Figure E] shows the evolution of the SFE for runs Mil-Mb.ll-Ma2e4, Mil-Mb.27-Ma2e4, 
Mil-Mb.02-Ma2e4 and Mi0-Mb.l0-Ma2e4. The first three runs differ only in the strength of 
the initial turbulence, measured by M SiDg d (cf. Tabled]). The last run has nearly the same 
value of M Sj b g d as Mil-Mb.ll-Ma2e4, but with no inflow velocity, in order to assess the 
amount of star formation induced solely by the turbulent field, in the absence of colliding 
streams. 

From these figures, we see that tsF becomes longer and (SFE) becomes smaller as larger 
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t[Myr] t[Myr] 

Figure 5. Evolution of the "absolute" SFE (SFE a t, s ,25i left panel) and the "relative" SFE (right panel) for four runs, character- 
ized by various values of the initial background turbulent Mach number M s ,b g di indicated by the entry "Mb.##" in the runs' 
names. Three of the runs have the same inflow Mach number, M\ n i = 1.25, while the fourth run has no inflows (A^i n f = 0), in 
order to assess the SFE due exclusively to the initial background turbulence. 



values of M s> b g d are considered. This appears to be due to the stronger fragmentation in- 
duced by the turbulent velocity field, which in the extreme case of Mil-Mb.27-Ma2e4 almost 
obliterates the inflows, and produces scattered clumps througout the simulation box, with 
very little remaining of the cloud formed by the inflows, as illustrated in Fig. El left panel. 

It is also worth noting that run Mi0-Mb.l0-Ma2e4, which has no inflows, has a much 
larger t$F and a lower (SFE) than run Mil-Mb.ll-Ma2e4, which differs from the former 
only in the presence of the inflows. As shown in Fig. El this run also produces scattered 
clumps throughout the numerical box, although not as profusely as Mil-Mb.27-Ma2e4. So, 
we conclude that sink formation is still dominated by the colliding streams in Mil-Mb.ll- 
Ma2e4, although a small fraction of the sinks is contributed by the global turbulence. 

Figure [7] summarises the results of this section. A clear trend of a decreasing SFE (seen 
in both SFE abs 2 5 and SFE rel)5 ) with increasing M s bgd is seen, which we interpret as a result 
of the reduction of the fragment mass with increasing turbulence strength (at constant total 



Ballesteros-Paredes et al. 



20061 ). and of the fact that smaller-mass fragments tend to 



have smaller SFEs ( §3.31) . 



3.3 SFE vs. inflow mass 

The last dependence of the SFE we analyse is on the mass content of the colliding inflows. 
So, we consider inflows of various radii. However, since a very narrow inflow is necessarily 
more poorly resolved, we consider smaller simulation boxes in two of the cases, in order to 
better resolve the resulting clouds. Specifically, as shown in Table [[J i?i n f in run Mil-Mb.02- 
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Figure 6. Face-on images (in projection) of Mil-Mb.27-Ma2e4 {left panel) and Mi0-Mb.l0-Ma2e4 (right panel) at the time each 
one is beginning to form stars. Note the scattered structure, due to the turbulence forming clumps throughout the numerical 
box. In run Mil-Mb.27-Ma2e4, in which the colliding streams are present, the turbulence almost completely obliterates the 
"main" cloud formed by the streams. 



Ma5e3 is half that in Mil-Mb.02-Ma2e4. Run Mil-Mb.06-Ma2e3 has R mi equal to that of 
Mil-Mb.02-Ma5e3, but half the length, since the numerical box size of the former is half that 
of the latter. Finally, Mil-Mb.06-Ma6e2 has R in{ equal to half that of Mil-Mb.06-Ma2e3. 
So, the total mass contained in the inflows of runs Mil-Mb.02-Ma2e4, Mil-Mb.02-Ma5e3, 
Mil-Mb.06-Ma2e3 and Mil-Mb.06-Ma6e2 is, respectively, 2.26 x 10 4 , 5.64 x 10 3 , 2.42 x 10 3 , 
and 6.04 x 1O 2 M . 

Figure [S] shows the evolution of SFE abs 2 5 and SFE rel 5 for these runs. We see that there 
is a general trend for the SFE, in both its forms, to increase with the total mass involved in 
the stream collision. There is only a reversal to this trend in the relative SFE between runs 
Mil-Mb.02-Ma2e4 and Mil-Mb.02-Ma5e3 because the latter has a large early maximum of 
SFE reli5 , although later it decreases, in a period of mass accumulation in the cloud at low 
SFR. Other than that, the trend is general, as shown in Fig. [9l 

It is important to note that for all runs in this series we used the same energy injec- 
tion rate of the turbulence driver. However, runs Mil-Mb.06-Ma6e2 and Mil-Mb.06-Ma2e3, 
performed in a smaller computational box, have an initial, background rms turbulent Mach 
number that is roughly 2.5 times larger than that of runs Mil-Mb.02-Ma2e4 and Mil-Mb.02- 
Ma5e3. This may additionally reduce the SFE because of the additional fragmentation it 
produces, but we see that the trend of the SFE to decrease with decreasing inflow mass 
holds generally even at the same physical box size, so the result appea rs robust. 

The trend discused above can be put in the context of the theory of iKrumholz &: McKee 
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Figure 7. Dependence of SFE a b Sj 25 an d SFE ro i 5 on the rms Mach number, M s ,bgdi OI the initial (background) turbulent 
velocity perturbations. Both indicators arc seen to decrease with increasing M s bgd as a consequence of the progressively 
stronger fragmentation induced by the turbulence. 
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Figure 8. Evolution of the "absolute" SFE (SFE abSj25 , left panel) and the "relative" SFE (right panel) for runs Mil-Mb.02- 
Ma2e4, Mil-Mb.02-Ma5e3, Mil-Mb.06-Ma2e3, and Mil-Mb.06-Ma6e2. The radius of the cylindrical inflows for these runs is 
respectively 32, 16, 16, and 8 pc. Runs Mil-Mb.02-Ma2e4 and Mil-Mb.02-Ma5e3 are performed in a 256-pc box, with inflow 
length 112 pc, while runs Mil-Mb.06-Ma6e2 and Mil-Mb.06-Ma2c3 arc performed in a 128-pc box with the same number of 
SPH particles (thus being better resolved) and a 48-pc inflow length. 
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mass is observed. 



( 120051 . hereafter KM05) for the SFRfj. This theory predicts a dependence of the SFRg on 
the virial parameter a and the rms Mach number M s . Here, 



a = 2Evn/\E, 



kin / I -^grav | > 



(2) 



where -Ekin = M.Av 2 /2 is the cloud's kinetic energy, M. is the cloud's mass, Av is its rms 
turbulent velocity dispersion, and E grav is the cloud's gravitational energy 

For our flattened clouds, we compute the gravitational energy assuming they can be 
approximated as infinitely thin, uniform disks of radius R, and write 



-Sgrav = / d 2 x = 27rS / rcj)(r)dr, 

J A JO 



(3) 



where A is the area of the disk, £ is the (unifor m) surface density, and cj) is the gravitational 



poten tial. In our case, the latter is given by (jWyse fc Mayall 
20M) 



1942 



Burkert fc Hartmann 



(r) = -AGT,RE(r/R), 



(4) 
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where E is the second complete elliptic integral. Thus, the gravitational energy is 

£ = -8nGJ: 2 R [ R rE(r/R)dr = -8nGJ: 2 R 3 f 1 xE(x)dx 

Jo Jo 

To compute the kinetic energy of the clouds, we note that the relevan t velocity dispersion 
is the one produced in the clou ds as a consequence of the inflow collision (IHeitsch et al.ll2005 ; 



Vazquez-Semadeni et al 



20061 ). rather than the initial background turbulent velocity, which 



is much smaller. Since all four simulations analysed in this section have the same v in f, we 
use the value of t> rms measured for Run Mil-Mb.02-Ma2e4 in § 13.11 (Fig. S]) for all of them, 
namely t> rms = 0.5 km s _1 . Noting that this value is a one-dimensional velocity dispersion, 
we take Av = \/3 v rms . 

We finally obtain, from equations (jSJ) and ((Sj), 
224nRAv 2 

a= 45GM • (6) 
Figure [TU] shows the results of this exercise. The solid line shows the simulation data, while 
the s traight dotted line sh o ws a least squares fit to them. The dashed line shows the result 



Krumholz fc McKed (120051 ). given by 



from_ 

where we have taken M s = Av/c s . We see that the prediction by KM05, although being 
numerically within the same range as the data, exhibits a significantly shallower slope. 
Specifically, the fit to our data has a slope —1.12 ± 0.37, where the uncertainty is the la 
error of the fit, while the slope of the KM05 prediction, —0.68, lies beyond this error. We 
discuss this result further in 



4 SUMMARY AND DISCUSSION 

In this paper, we have considered the scenario of molecular cloud formation by WNM stream 
collisions, and investigated the dependence of the SFE on three parameters of this scenario, 
namely the inflow speed, the rms Mach number of the background medium, and the total 
mass contained in the inflows. Since the SFE, defined as in eq. (pQ), is a time- dependent 
function because the cloud continues to accrete mass from the WNM while it forms stars, 
we have considered two estimators of its time integral, namely the absolute SFE after 25 
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Figure 10. Dependence of SFE rc i 5 on the virial parameter a. The straight dotted line shows a le ast squares fit to the simulation 
data, with slope — f.12, while the dashed line shows the result from [Krumholz &: McKeel {2005), with slope —0.68. 

Myr from the start of the simulation, SFE aDS ,25> and the "relative" SFE, 5 Myr after the 
onset of SF in the cloud, denoted SFE re i i5 . 

We have found a wide range of values of these estimators as we vary the parameters of 
the simulations. In general, the SFE decreases, although moderately, with increasing inflow 
velocity. In particular, SFE re i 5 decreases from ~ 0.4 to ~ 0.2 as M s in f increases from 1.25 
to 3.5. Note that the runs in this case have similar SFRs, and the decrease in the SFE and 
(SFE) is due to the faster increase in cloud mass rather than to a decrease in the SFR. That 
the SFR is similar in all three runs, in spite of the larger gas mass is probably due to the 
larger turbulent velocity dispersion in the dense gas caused by the larger inflow speed, which 
tends to inhibit the SFR, thus compensating the tendency to have a larger SFR due to the 
larger cloud masses. 

Similarly, the SFE decreases with increasing background turbulence strength, as the 
latter progressively takes a dominant role in the production of the dense gas but, due to the 
relatively small scales at which the turbulence is excited, the clouds and clumps formed by 



16 Rosas- Guevara et al. 

it are significantly smaller than the cloud formed by the coherent stream collision. In this 
case, SFE re i,5 decreases from ~ 0.4 to ~ 0.03 as the rms Mach number of the background 
turbulence increases from ~ 0.02 to ~ 0.3. 

Finally, the SFE in general decreases with decreasing mass of the inflows (at constant 
^inf). This may be a consequence of the fact that clouds formed by the collision of our inflows 
always have roughly the same density, temperature, and velocity dispersion, so s maller clouds 



are m ore weakly gravitationally bound, a condition known to decrease the SFE ( jClark et al 



20051 ). The end result is that SFE re i 5 decreases from ~ 0.4 to ~ 0.03 as the mass in the 
colliding streams decreases from ~ 2.3 x 1O 4 M to ~ 6OOM . It is important to stress that 
this result is not at odds with the well known fact that the SFE increases as the object mass 
decreases from the mass scale of a giant molecular cloud (Ai ~ 10 4 - 1O 6 M , SFE ~ 0.02; 
Myers et al. 1986) to that of a cluster-forming core (M ~ 1O 3 M ; SFE ~ 0.3 - 0.5, Lada & 
Lada 2003), because in this case the cores' mean densities are much larger than those of the 
GMCs, while in our case the mean densities of the various clouds are always comparable. 
Thus, our clouds do not conform to Larson's (1981) density-size scaling. 

The latter results, expressed in terms of the virial parameter a, exhibit partial agree- 
ment with the prediction by KM05 for the dependence of the SFE after a free-fall time 
(what those authors called the star formation rate per free-fall time, or SFRg), which is 
directly comparable to our SFE re i 5 . Although their prediction, without any rescaling, lies 
in the same range of values as our observed efficiencies, it contains a much shallower de- 
pendence on a than we observe. This may be due to the fact that those authors assumed 
that the clouds were supported by turbulent pressure, possibly provided by stellar feedback, 
while our simulations l ack such support. However, the notion of turbulent support has been 



questioned recently by I Vazquez- Semadeni et al.l (120081 ). and the low observed efficiency may 



be th e result of the dispersal (rather than support) of the parent cloud by its stellar prod 



ipp 
irJL 



ucts ( lHartmann. Ballesteros-Paredes fc Berginll200ll ). More work is needed to determine the 



causes of the discrepancy. However, it is noteworthy that Run Mil-Mb.06-Ma2e3, which has 
a value of a closest to unity, has a reasonably realistic value of SFE re i 5 ~ 4%. 

We conclude that the SFE, even in the absence of further agents such as magnetic 
fields and stellar feedback, is a highly sensitive function of the parameters of the cloud 
formation process, and may be responsible for significant intrinsic scatter in observational 
determinations of the SFE. 
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Table 1. Run parameters 



Run 


Run 


^box 


^ i nf 


L mi 




-"^s, bga 


1 L int 


* v 1 box 


J v 4111 


J v >■ part 


number 


name 


[PC] 


[PC] 


[ km s- 1 ] 






[PC] 


[M e ] 


[M Q ] 


[M@] 


1 


Mil-Mb.ll-Ma2c4 


256 


112 


9.20 


1.25 


0.11 


32 


5.25 x 10 5 


2.26 x 10 4 


0.32 


2 


Mi2-Mb.ll-Ma2e4 


256 


112 


18.41 


2.50 


0.11 


32 


5.25 x 10 5 


2.26 x 10 4 


0.32 


3 


Mil-Mb.27-Ma2e4 


256 


112 


9.20 


1.25 


0.27 


32 


5.25 x 10 5 


2.26 x 10 4 


0.32 


4 


Mil-Mb.02-Ma2e4 


256 


112 


9.20 


1.25 


0.021 


32 


5.25 x 10 5 


2.26 x 10 4 


0.32 


5 


Mi2-Mb.02-Ma2e4 


256 


112 


18.41 


2.50 


0.024 


32 


5.25 x 10 5 


2.26 x 10 4 


0.32 


6 


Mi3-Mb.02-Ma2e4 


256 


112 


25.77 


3.50 


0.025 


32 


5.25 x 10 5 


2.26 x 10 4 


0.32 


8 


Mil-Mb.02-Ma5c3 


256 


112 


9.20 


1.25 


0.020 


16 


5.25 x 10 5 


5.64 x 10 3 


0.32 


10 


Mil-Mb.06-Ma6e2 


128 


48 


9.20 


1.25 


0.057 


8 


6.57 x 10 4 


6.04 x 10 2 


0.04 


11 


Mi0-Mb.l0-Ma2e4 


256 


112 


0. 


0. 


0.10 


32 


5.25 x 10 5 


2.26 x 10 4 


0.32 


12 


Mil-Mb.06-Ma2e3 


128 


48 


9.20 


1.25 


0.058 


16 


6.57 x 10 4 


2.42 x 10 3 


0.04 



